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SUMMARY 

Flight  data  from  a 60°  delta  wing  aircraft  have  been  analysed  using  a modified 
Newton-Raphson  parameter  estimation  procedure.  The  model  equations  used  for  the 
analysis  were  extended  to  account  for  sideslip  vane  errors  and  for  lateral  accelerometer 
position  error.  Lateral  derivatives  extracted  from  the  data  have  been  compared  with  wind 
tunnel  measurements  and  theoretical  estimates  and  areas  of  agreement  and  disagreement 
identified.  The  method  has  also  been  applied  to  the  analysis  of  fin  loads  measured  in  flight 
and  some  tentative  conclusions  reached.  The  results  confirm  the  effectiveness  of  the  para- 
meter identification  procedure  in  flight  test  analysis  and  its  ready  applicability  to  a variety 
of  related  problems.^ 


POSTAL  ADDRESS:  Chief  Superintendent,  Aeronautical  Reaearch  Laboratoriea, 
Box  4331 , P.O.,  Melbourne,  Victoria,  3001 , Australia. 
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1.  INTRODUCTION 


The  modified  Newton-Raphson  parameter  estimation  method  described  in  References  1 
and  2 has  been  adapted  for  use  on  the  ARL  PDPIO  computer.  In  support  of  mathematical 
modelling  activities  being  presently  undertaken  the  method  provides  a means  for  systematic  and 
rational  validation  of  a mathematical  model  as  well  as  a means  of  obtaining  aerodynamic  data 
not  otherwise  available.  In  Reference  3 a detailed  study  of  the  method  was  made  to  establish 
requirements,  such  as  input  design,  sampling  rates,  record  length,  etc.,  for  a successful  appli- 
cation of  the  method  to  the  analysis  of  response  records.  The  method  was  subsequently  applied^ 
to  the  analysis  of  actual  flight  test  data  of  the  longitudinal  response  of  a delta  wing  fighter  aircraft 
to  elevator  inputs.  The  present  study,  looking  at  the  lateral  response  of  the  aircraft  to  rudder 
inputs,  raises  a number  of  new  aspects  which  complement  the  work  of  Reference  4.  In  addition 
the  analysis  is  extended  to  the  estimation  of  airload  parameters  which  govern  the  forces  experi- 
enced by  the  fin  when  the  aircraft  is  in  flight. 

The  basic  concept  of  the  estimation  method  is  summarised  in  Figure  1 . The  objective  is  to 
determine  that  value  of  the  mathematical  model  parameter  vector,  c,  which  minimises  a cost 
functional,  J,  proportional  to  the  difference  between  the  measured  response,  z,  and  the  calculated 
response,  y,  i.e. 

J = (1/N)  Z (z<  - y<)^Z)i(z<  - y<)  -I-  (c  — — Co)  (1) 

/=i 

where  Di  and  D2  are  diagonal  weighting  matrices  and  the  second  term  is  a weighted  mean  square 
difference  between  the  parameter  vector  c and  a stipulated  a priori  value,  Co. 

The  mathematical  model  is  described  by  a state  equation  of  the  form 

X = -h  Ru  (2) 

and  the  calculated  response  is 

y = Fx  -I-  Cu  + b.  (3) 

The  measured  response  is  assumed  to  be  given  by 

2 = y -|-  n f4) 

where 

X is  the  vector  of  state  variables, 
u is  the  vector  of  control  variables, 
y is  the  response  vector, 
b is  a measurement  bias  vector, 

2 is  the  measured  response  vector, 
u is  a random  noise  vector, 

A,  B,  Fand  G are  matrices  of  parameters  defining  the  model. 

The  parameter  vector,  c,  contains  some  or  all  of  the  elements  of  A,  B,  Fand  G,  the  elements 
of  b and  the  initial  conditions.  The  iterative  modified  Newton-Raphson  procedure  is  used  to  find 
that  value  of  c which  minimises  the  cost,  J.  An  indication  of  the  quality  of  the  estimates,  c,  is 
given  by  the  Cramer-Rao  bound,  acn,  which  is  a lower  bound  on  the  variance  of  the  estimates 
(see  Ref.  2). 

The  detailed  form  of  the  model  equations  is  developed  in  Section  2 and  includes  allowance  for 
sideslip  vane  errors  and  lateral  accelerometer  corrections.  In  addition,  equations  for  the  fin  loads 


I 


are  also  outlined.  Section  3 describes  the  flight  data,  pre-analysis  data  conditioning,  and  discusses 
the  a priori  values  and  weights  used  in  the  analysis.  Results  are  presented  in  Section  4,  together 
with  discussion  and  comparison  with  theoretical  values  of  several  of  the  aerodynamic  derivatives. 


2.  MATHEMATICAL  MODEL 

In  this  section  the  detailed  form  of  the  mathematical  model,  described  in  general  by 
Equations  (I)  and  (2),  is  developed.  Firstly,  the  basic  linear  small  perturbation  equations  for 
aircraft  lateral  response  to  control  input  are  outlined.  Next  an  extension  of  the  model  is  developed 
in  order  to  deal  with  errors  in  sideslip  vane  measurements  and  this  is  followed  by  a method  for 
incorporating  lateral  accelerometer  corrections  into  the  model  for  a case  where  these  corrections 
cannot  be  directly  made  to  the  measurements.  Finally  a model  for  obtaining  fin  load  parameters 
is  described. 

2.1  Basic  Model 

The  aircraft  is  assumed  to  be  initially  trimmed  in  steady  level  flight  at  speed  Ve,  incidence 
cte  and  attitude  Se-  Lateral  small  disturbance  motions  about  this  trim  state  can  be  produced  by 
either  aileron  or  rudder  inputs.  The  linearised  equations  describing  the  small  disturbance  lateral 
motions  following  a control  input  can  be  written,  in  body  axes,  as  follows  (see  Reference  S for 


details): 

mv  = YfP  + Y,,Br  + Yi^&a  — muer  + mwep  + mg  cos  ^ (5a) 

Ixp  — Ixzf  = LfP  + Lpp  -1-  Lrr  + L0P  Lf^Sr  -t-  (5b) 

Itf  — Ixtp  = + Npp  -\-  Nrr  + -|-  Nffir  Nt^Bo  (5c) 

<i>  = P (5d) 

where  w*  = sin  ««  ~ K«ae  and  m*  = Ve  cos  ««  Ve  for  small 

Further,  from  the  definition  v = 1^  sin  it  follows,  for  constant  speed  and  small  sideslip,  /3, 

that 

V = Vep.  (6) 


Hence  Equation  (5a)  becomes  an  equation  for  /3  by  substitution  of  Equation  (6).  This  can 
then  be  used  to  eliminate  the  jS  terms  from  Equations  (5b)  and  (5c).  In  the  present  case,  however, 
this  is  not  necessary  since  it  can  be  assumed  that  Lf  and  Nf  are  negligible.  After  some  algebra. 
Equations  (5b)  and  (5c)  can  be  rewritten  as 

p = L' ^ -f-  L' pp  + L'rr  + L'lf&r  -f-  L'l^Sa  (7a) 

r = N’fp  + N'pP  + N'rr  + N',J&r  -|-  N',M  (7b) 

where  the  dashed  derivatives  are  defined  as 

L'c  = I'xzNe  -I-  LcH'x  (8) 

N'e  = I'xeLe  + Nc/l'e 

with  the  subscript  c standing  for  p,  r,  8r  or  8a  and 

I'x  = Uxlz  - Ihz)//z 

I'z  = {Ixlz  - PxzMIx  (9) 

I'xz  = txzKIxlz  — f^xz)- 

Since  /«  is  small  in  the  present  case  and  can  be  assumed  to  be  zero  the  relations  given  by 
Equations  (8)  and  (9)  are  considerably  simplified,  e.g.  L'e  = LdU,  N'c  = Ntih- 
Finally,  the  equation  for  the  lateral  acceleration,  rip,  follows  from  the  definition 

iglVe)ny  = (u  + Uer  — HVp  — g COS  de-  (10) 

~ iYflmVe)^  + {YJmVe)8r  + {Y,JmVe)8a. 
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Thus  the  basic  mathematical  model  used  for  the  analysis  of  the  lateral  motion  may  be 
summarised  in  matrix  form  as  follows. 


The  state  equations  are; 


p 

~ Yf 

mVo 

-1 

g cos  6,  ~ 

y* 

p 

" y,r 

mV, 

Yta 

mV, 

Yo 

Sr 

p 

= 

hi 

hi 

Ix 

Lr 

I, 

0 

p 

+ 

L,, 

Ix 

L,. 

Ix 

Lo 

s. 

Np 

Nr 

N>r 

No 

1 

It 

h 

h 

u 

r 

It 

h 

J. 

0 

1 

0 

0 _ 

0 

0 

0 _ 

The  equations  for  the  response  variables  are: 
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(12) 


The  parameter  vector  is  given  by: 

c = [YfImVi,  YtJmVt,  Y^JmVf,  Lgjlx,  Lpjlx,  I-rlIxt  l-Srl^x,  LtJIx, 

Nflh,  Nplh,  Nrih,  Ntrlh,  NJh,  Yo,  U No,\^  (13) 

The  measurement  bias,  b (appearing  in  Equation  (3)),  has  been  neglected  in  the  above 
equations  [and  the  parameters  To,  Lo,  and  No  have  been  introduced  to  account  for  possible 
state  bias. 

• ■ 

2.2  Correctkms  to  Sideslip  Vaae  Measarements’ 

The  sideslip  measurements,  jS*,,  are  influenced  by  yaw  and  roll  rate  effects,  and  are  related 
to  the  true  sideslip,  by  the  following  equations  (Reference  5): 

= (14) 

with 

= Pi^vlYt)  - r(xv/Ve) 

where  Xv  and  Zv  are  the  distances  of  the  vane  forward  and  below  the  centre  of  gravity.  In  order  to 
account  for  possible  calibration  slope  errors.  Equation  (14)  has  been  modified  to 

p„  = GP-p,  (15) 

with  G becoming  a parameter  to  be  identified  during  the  matching  process. 

Since  the  vane  was  not  dynamically  balanced,  it  also  responds  to  accelerations.*  The  cor- 
rection to  the  measured  sideslip  due  to  lateral  acceleration  at  the  vane  can  be  written: 

Ps  = mXolVe)  - piZvIVe)  + n/gIVo))  (16) 

where  yaw  acceleration,  f,  and  roll  acceleration,  p,  are  in  radians/sec^  and  lateral  acceleration, 
np,  is  in  g’s.  An  estimated  value  for  the  constant  of  proportionality,  K,  depending  on  vane  charac- 
teristics, is  given  in  Section  3.1.  The  measured  sideslip.  Pm,  can  now  be  related  to  the  true  side- 
slip, p,  by  the  equation 
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¥ 

ii 

Pm  = Gp-p.+Pa  (17) 

or 

Pm+P.  = GP  + mx^jVe)  - piz^lVe)  + n^/Ve)).  (18) 

By  substituting  for  r,  p,  and  n»  from  Equations  (7)  and  (10),  and  assuming  8*  to  be  zero. 
Equation  (18)  can  be  expanded  to  give: 

Pm+P.  = [G+  K[(Xviye)Nflh  - {z^lV,)LflU  + Yf/mVeW 

+ K[{Xviyt)Nvlli  ~ i.Zviyt)LplIx\p 

+ K[{Xviy.)Nrlh  - (Zviye)Lrllx]r 

+ Kl(Xviye)N,JU  - iZx/ye)LM&r 

= fli/S  + a2p  + Oar  + biBr  (19) 

The  corrections  due  to  roll  and  yaw  rates,  /3,,  are  readily  calculated  from  Equation  (14)  since 
measurements  of  p and  r are  available.  Thus  (/Sm  + ^»)  rather  than  can  be  treated  as  the 
response  variable  to  be  matched  with  the  model  resulU  given  by  the  right  hand  side  of 
Equation  (19). 

23  CorrectloHS  to  Lateral  Accelerometer  Measurements 

Since  the  lateral  accelerometer  is  not  located  at  the  centre  of  gravity  the  acceleration 
measured  by  it  will  differ  from  the  value  at  the  centre  of  gravity.  If  the  accelerometer  position  is 
given  by  the  body  axes  co-ordinates  (xo,  ya,  za)  then  the  relation  between  the  measured  and  c.g. 
values  can  be  written^: 

”ym  = »v  + (w  — ^«P  - >'»'■*  — (20) 

where  ny„  is  the  measured  acceleration  and  riy  is  the  value  at  the  c.g.  The  dominant  correction 
term  in  the  present  investigation  was  found  to  be  that  due  to  roll  acceleration,  p,  caused  by  the 
location  of  the  lateral  accelerometer  approximately  0 - 6 m (2  ft)  below  the  c.g.  Hence  Equation  (20) 
can  be  simplified  to 

nym  = Oy  - (Zalg)p.  (21) 

If  measurements  of  p were  available  then  the  correction  (zalg)p  could  immediately  be  applied 
to  the  measured  acceleration.  However,  as  p is  not  available  in  the  present  case  one  possibility 
is  to  proceed  with  the  identification  without  correcting  and,  using  the  value  of  p thus  calculated, 
then  make  the  correction  and  repeat  the  matching  process.  Another  alternative,  adopted  here,  is 
to  allow  for  the  correction  within  the  mathematical  model  itself  by  calculating  the  right  hand  side 
of  Equation  (21)  and  matching  this  directly  with  the  measured  value,  ny„.  Substituting  for  ny 
and  p in  Equation  (21)  from  Equations  (7)  and  (10),  and  assuming  8*  to  be  zero,  results  in 

(giye>y„  = [Yf/mye  - (Zaiye)LflIx]P  - {ZaiyeY^LplIx)? 

— (Zo/KeKWAt)''  + [Ytf/mye  — iZaiye)LiJIx]br 

= iyfilr»ye)P  + {yplmye)p  + (yr/mK,)r  + (y«,/mK*)Sr  (22) 

where  y^,  yp,  yr  and  yt,  are  modified  side  force  derivatives. 

Thus  Equation  (22),  for  ny„,  replaces  Equation  (10),  for  riy,  in  the  mathematical  model.  The 
full  equations,  including  the  corrections  outlined  in  this  section  and  the  sideslip  corrections  of 
Section  2.2  may  now  be  summarised.  Since  there  was  no  aileron  input  in  the  course  of  the  flight 
tests  the  equations  have  been  simplified  by  removing  the  aileron,  8o,  term  from  the  input  vector. 
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The  state  equations  are: 
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The  response  equations  become; 
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(24) 


where  bo,  yo  are  introduced  to  account  for  possible  bias.  The  parameter  vector,  excluding  bu 
ypImVt  and  yrImVe,  which  were  held  constant  throughout  as  discussed  in  Section  3.3,  can  be 
written 


c = lY,lmVe,  Y,JmVe,  y^ImVe,  y,,lmV„  L^jh,  LplU,  Lrih,  LJU, 

bfp/Ji,  Nrjlz,  Nffjlz,  fll,  Oi,  <73,  Yo,  Lo,  No,  bo,  Yo]^ 


(25) 


2.4  Equations  for  Fin  Loads 


While  the  response  vector.  Equation  (24)  in  Section  2.3,  consists  of  elements  usually  measured 
in  the  course  of  stability  and  control  flight  testing,  the  method  of  analysis  does  not  place  any 
restrictions  on  the  elements  making  up  the  response  vector.  For  example.  Reference  7 applies 
the  method  to  the  determination  of  in-flight  airload  parameters.  During  the  flight  tests  under 
consideration  here,  strain  gauges  in  the  fln  measured  Shear,  S,  and  Tension,  T,  forces.  Following 
Reference  7,  these  forces  can  be  related  to  the  state  and  control  variables  through  the  following 
linearised  small  disturbance  equations: 

S = SpP  -f-  Spp  -f-  Sfr  -f"  St,Sr  (26) 

T=Tffi  + Tpp  + Trr  + Tt,hr  (27) 

In  Equations  (26)  and  (27)  the  inertial  and/or  gravitational  contributions  to  S and  T have 
been  neglected  as  being  small.  To  identify  the  Shear  and  Tension  parameters  (S^,  Tp,  etc.)  the 
response  vector  can  be  readily  extended  to  include  S and  T as  additional  elements  so  that  the 
S and  T records  may  be  matched  simultaneously  with  the  other  records.  Alternatively,  a two 
part  procedure  can  be  followed  since  the  airload  parameters  of  Equations  (26)  and  (27)  are 
completely  uncoupled  from  the  other  aircraft  parameters.  The  first  step  would  identify  the  air- 
craft parameters  as  usual  and  these  would  subsequently  be  held  constant.  In  the  second  stage 
only  5 and  T would  be  matched.  The  mathematical  model  for  the  second  stage  would  have  the 
same  state  equations  as  given  by  Equation  (23)  while  the  response  equations  would  be; 
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where  So,  To  are  bias  terms.  The  parameter  vector  becomes 

c — Sp,  Sr,  Sf^,  Tfi,  Tp,  Tr,  Tf„  So,  T’o]*'  (29) 

The  second  method  has  been  used  in  obtaining  the  results  described  in  Section  4.3.  No 
significant  change  in  the  accuracy  of  the  parameters  is  expected  compared  to  the  method  which 
treats  all  equations  simultaneously,  but  a significant  reduction  in  computer  time  results  (see 
Reference  7). 

3.  FUGHT  DATA 

The  flight  data  to  be  analysed  is  discussed  in  this  section.  The  test  conditions  and  the  form 
of  the  available  data  is  first  outlined  followed  by  a description  of  the  corrections  and  pre- 
processing applied  to  the  basic  data  in  preparation  for  the  Modified  Newton-Raphson  program. 
Finally  the  a priori  values  used  for  each  of  the  parameters  of  interest  is  summarised  followed  by 
a discussion  of  the  values  chosen  for  the  elements  of  the  weighting  matrices  Di  and  Dz. 


3.1  Test  Conditions 

The  flight  test  data  consisted  of  time  response  records  of  a 60°  delta  wing  fighter  aircraft 
following  a pilot  applied  rudder  pulse  input.  Records  for  two  different  Mach  numbers  were 
available,  namely  M = 0-96  and  M = 0 72.  The  relevant  test  conditions  for  these  flights, 
referred  to  as  flights  RKl  and  RK2  are  summarised  in  Table  1. 


TABLE  1 

FUgkt  Test  ConditioM 


Details 

• « 

. . .Flight  RKl 

Flight  RK2 

M 

0-96 

0-72 

h 

9076  m (29,770  ft) 

2985  m (9791  ft) 

Ve 

292  -4  m/s  (568  kt) 

236  m/s  (459  kt) 

3-73° 

2-71° 

Oe 

2-87° 

3-12° 

(8«)« 

-0-78° 

-0-61° 

c-g. 

50-33%  of  ^ chord 

48 -95%  of  4.  chord 

K 

0-293 

0-237 

The  values  of  K given  in  Table  1 were  estimated  using  the  formula 

K = 0- 00984  Volg  ^30) 

which  was  derived  in  Reference  4. 

Measured  data  consisted  of  response  time  histories  of  /3,  p,  r,  rip,  and  &r-  Approximately 
IS  seconds  of  record  was  available  but  only  S seconds  were  used  in  any  single  analysis.  This 
represented  approximately  2 to  3 cycles  of  oscillation.  The  data  rate  was  60  samples/sec  and  from 
inspection  of  the  flight  records  it  was  determined  that  the  natural  frequency,  cop,  was  about  3 • 2 
rad/s  for  flight  RKl  and  3'7  rad/s  for  flight  RK2.  Assuming  an  in-between  value  of  3 5 rad/s 
this  implies 


1 

1 
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a),A/  = 3-5/60  = 0-06 
u.„r  = 3-5x5=  17-5 


(31) 


where  Ar  is  the  sampling  interval  and  T the  total  record  length. 

In  Reference  3 it  was  found  that  values  of  unAt  <0-14  and  to^T  > 14  were  adequate  for 
matching  of  the  response  to  a rudder  doublet  input.  Although  in  the  present  case  the  input  is  a 
pulse  rather  than  a doublet,  the  values  given  in  Equation  (31)  suggest  that  the  sampling  rate  and 
record  length  should  be  ample  for  successful  matching.  In  fact  no  difficulties  of  convergence 
or  identification  were  encf>.ntered.  Normally  seven  iterations  were  specified  in  the  program, 
though  it  was  found  that  good  convergence  was  invariably  achieved  by  the  fourth  or  fifth  itera- 
tion with  only  small  improvements  subsequently. 

3.2  Data  Conditions 

Before  proceeding  with  the  analysis  a certain  amount  of  pre-processing  of  the  data  was 
performed.  This  included  the  following— 

(a)  Roll  and  yaw  rate  gyro  misalignment  corrections  according  to  the  equations  given  in 
Reference  5.  Using  the  gyro  alignment  data  the  corrected  roll  rate,  p,  and  yaw  rate,  r, 
are  obtained  from  the  measured  roll,  pitch  and  yaw  rates,  p„,  qm  and  r*  respectively, 
according  to  the  equations 

p = 0-9996/»„  -h  0-00131^m  - 0-0291r„ 
r = 0 ■ 0278/>™  - 0 - 04524r„  1 - 00092r«  (32) 

(b)  Corrections  to  angle  of  sideslip  measurements  due  to  yaw  and  roll  rate  effects  as  dis- 
cussed in  Section  2.2.  The  corrected  sideslip  angle.  Pm  + Pi,  was  subsequently  matched 
where  Pm  is  the  measured  value  and  p,  is  the  correction  given  by  Equation  (14). 

< (c)  A smooth  curve  was  hand  fitted  to  the  rudder  input  data  as  shown  in  Figures  2a  and  2b. 

Although  this  smoothing  is  not  strictly  necessary  it  was  thought  desirable  to  remove 
the  obvious  spikes  caused  by  a noisy  transducer  (see  Figure  2a  in  particular)  since  the 
identification  method  makes  no  allowance  for  noise  in  the  input.  For  a relatively  noise 
free  record,  such  as  that  in  Figure  2b,  the  effect  of  smoothing  on  the  results  is  hardly 
discernible. 

(d)  Because  the  roll  rate  gyro  saturated  at  a fairly  low  level,  the  initial  peak  of  the  roll  rate 
record  was  cut  off  as  shown  in  Figures  3a  and  36.  The  missing  portion  was  filled  in  by 
hand  after  choosing  an  apprqprjate  peal^  level  by  comparjsoa  with  the  ether  records.. 
The  result  is  also  shown  in  Figures  3a  and  36.  The  effect  of  this  rather  arbitrary  procedure 
on  the  matched  results  will  be  discussed  further  in  Section  4. 

3.3  A Priori  Values  and  Weighting  Matrices 

The  a priori  values  assigned  to  the  various  parameters  were  mainly  obtained  from  wind 
tunnel  tests  (e.g.  Ref.  8)  or  else  were  estimated  from  data  sheet  sources.  Table  2 summarises  the 
values  used  in  the  two  cases  analysed. 


7 


TABLE  2 

A Priori  Parameter  Values 
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I 
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Parameter 

Flight  RKl 

Flight  RK2 

YfImVe 

-0184 

-0-292 

YJmV, 

-0  026 

-0-040 

y0lmVt 

-0120 

-0-226 

ypImVe 

0 00403 

0-00465 

yrlmVe 

-0'O139 

-0-00187 

ytrImVe 

-0015 

-0-021 

hUx 

-30-67 

-25-82 

Lpjlx 

-1-92 

-1-78 

Lrjlx 

0-661 

0-684 

Ltrlfx 

-5-38 

-7-57 

N0IU 

12-06 

17-44 

Npllz 

-0-0066 

0-0126 

Nrih 

-0-594 

-0-563 

N,Jh 

4-007 

5-79 

ai 

1-0076 

1-079 

as 

0-0012 

0-0015 

as 

-0-0059 

-0-0074 

hi 

0-033 

0-051 

Yo 

0 

0 

yo 

0 

0 

Lo 

0 

0 

No 

0 

0 

bo 

0 

0 

The  yawing  moment  derivatives,  Nf  and  have  been  adjusted  to  the  correct  centre  of 
gravity  position  for  each  flight.  The  derivatives  yp,  yr,  and  y,,  are  defined  by  Equation  (22) 
while  the  parameters  ai,  02,  as,  and  hi  are  defined  by  Equation  (19).  In  both  cases  the  a priori 
values  were  calculated  from  the  respective  definitions  using  Xv  = 9 ■ 1 5 m (30  ft),  z,  = 0 - 6 m (2  ft) 
and  Za  = 0-6  m (2  ft).  For  flight  RKl  the  values  of  the  aerodynamic  derivatives  appearing  in  the 
definitions  were  set  at  their  respective  a priori  values,  as  in  Table  2.  For  flight  RK2  on  the  other 
hand,  an  initial  run  was  made  with  Za,  Zv  and  Xv  all  set  to  zero  and  the  aerodynamic  derivatives 
thus  identified  were  used  to  evaluate  the  formulae  of  Equations  (19)  and  (22).  Several  of  the 
parameters  listed  in  Table  2 are  weak  in  the  sense  that  they  only  slightly  influence  the  calculated 
results  and,  conversely,  the  information  available  in  the  flight  records  is  insufficient  to  alter  them 
significantly  from  the  a priori  values  assigned.  This  applies  in  particular  to  yp,  yr  and  hi,  which 
were  consequently  held  constant  at  their  a priori  values  and  also  to  as  and  as  which  were  only 
allowed  to  vary  in  one  case  (see  Section  4.1). 

The  weighting  assigned  to  the  a priori  values  of  each  of  the  parameters  which  were  allowed 
to  vary  was  established,  as  in  Reference  4,  by  assuming  a standard  deviation,  a,  equal  to  20% 
of  the  parameter  value.  Then  the  relevant  element  of  the  weighting  matrix  Ds  becomes  simply 
l/c2.  The  choice  of  20%  is  not  critical  with  only  small  differences  in  the  identified  results  following 
from  a value  of  10  %.  In  general  the  smaller  the  value  of  a the  more  the  identified  parameter  values, 
especially  those  of  the  weaker  parameters,  tend  towards  the  a priori  value.  Those  parameters  in 
Table  2 with  zero  a priori  values  were  assigned  very  small  weighting  thus  allowing  them  freedom 
to  take  up  whatever  value  they  may  tend  towards,  no  matter  how  weakly.  In  addition,  during 
the  analysis  of  fin  loads,  all  the  parameters  in  Equation  (29)  describing  the  fin  loads  were  assigned 
zero  a priori  values  and  effectively  zero  weighting. 

The  elements  of  Di  determine  the  weighting  to  be  assigned  to  each  of  the  time  histories 
being  matched.  In  the  present  case  these  include,  at  one  time  or  another,  Pm  + Pt,  P,  r,  {.glVt)ny 
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as  well  as  shear,  S,  and  tension,  T,  records.  In  order  to  give  approximately  equal  weighting  to 
each  matched  time  history  a standard  deviation,  am,  equal  to  5%  of  the  maximum  value  of  each 
measured  record  was  selected  and  the  relevant  element  of  Di  set  to  l/o^m-  The  resulting  weights 
are  summarised  in  Table  3 below. 


TABLE  3 

ElemcBb'  of  Wcigbtiag  Matrix  D| 


Response  vector  element 

Flight  RKl 

night  RK2 

P 

0 001 

0 001 

P 

816 

3906  (657) 

r 

6400 

3906 

0 001 

0 001 

(.Pm  + P*) 

51653 

41649 

(giye)"y 

2mm 

1041233 

s 

I7I 

140 

T 

5-2 

3-8 

It  should  be  noted  that  the  final  identified  results  do  not  depend  critically  on  the  weights. 
Thus  for  flight  RK2  the  same  weight  (3906)  was  assigned  to  both  p and  r records  (Table  3) 
instead  of  assigning  to  p the  more  appropriate,  according  to  the  previous  argument,  value  of 
657.  The  only  discernible  difference  thus  produced  was  in  the  weighted  mean  square  match 
error  5®(p),  as  is  to  be  expected.  For  records  not  being  matched  the  corresponding  element  of  Di 
can  be  set  to  a small  value.  This  applies  to  p and  <f>  in  Table  3.  In  the  case  of  <f>  no  measurements 
were  available  for  matching  while  in  the  case  of  p the  corrected  value,  Pm  + Pt,  rather  than  P 
itself,  is  being  matched. 

4.  RESULTS 

In  this  section  the  results  of  the  parameter  estimation  analysis  of  the  flight  manoeuvres  are 
presented.  The  matches  obtained  between  the  mathematical  model  calculations  and  the  flight 
records  are  first  discussed.  This  is  followed  by  a closer  look  at  the  values  obtained  for  the  various 
aerodynamic  derivatives  and  a comparison  with  the  a priori  estimates  and  theoretical  estimates. 
Finally,  the  results  of  the  fin  load  analysis  are  presented. 

4.1  Results  of  Matching 

The  results  from  flight  RKl  at  A/  = 0-96  are  presented  in  Figures  4-7  and  in  Table  4 below. 
Table  4 summarises  the  a priori  values  and  respective  weights  in  the  first  column  with  the  results 
(plus  Cramer-Rao  bounds)  of  several  identification  runs  in  subsequent  columns.  Note  that  the 
parameters  az  and  <73,  though  usually  fixed,  were  allowed  to  vary  in  case  RKI-4.  As  can  be  seen 
from  Table  4,  the  identified  values  of  az  and  <73  did  not  differ  from  their  respective  a priori 
estimates  nor  did  the  Cramer-Rao  bounds  differ  from  the  a priori  weights.  It  is  clear  then  that 
<72  and  <73  are  weak  parameters  and  holding  them  fixed  in  all  other  cases  is  not  a restriction. 

The  rudder  input  in  each  case  is  the  smoothed  curve  shown  in  Figure  2a.  It  should  be  pointed 
out  that  smoothing  of  the  rudder  input  in  the  present  case  makes  little  difference  to  the  results. 
In  addition,  the  slight  difference  in  the  smoothed  and  measured  record  for  / > 1 • 30  sec  in  Figure 
2<i  is  not  significant.  Cases  have  been  run  with  8r  set  to  0004  rad.  for  t > I - 30  sec  with  almost 
identical  results. 

In  case  RKI-1,  shown  in  Figure  4,  all  the  available  flight  records  have  been  matched  using 
the  full  mathematical  model  summarised  in  Section  2.3.  The  match  is  not  as  good  as  may  be  hoped 
for  with  fairly  large  mean  square  errors,  i*,  as  shown  in  column  RKI-I  of  Table  4.  One  possible 
cause  for  this  may  be  inadvertent  errors  introduced  by  extrapolation  of  the  p record  as  shown 
in  Figure  3a.  To  check  this  possibility  the  match  was  repeated  with  the  p record  excluded.  The 
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results  arc  shown  in  Figure  5 and  column  RKI-2  of  Table  4.  Although  the  fits  to  the  sideslip, 
yaw  rate,  and  lateral  acceleration  records  have  improved  somewhat  as  shown  by  the  reduced 
mean  square  errors  (Table  4),  they  are  still  not  good.  At  the  same  time  the  removal  of  the  p 
record  has  removed  much  of  the  rolling  moment  information  so  that  the  rolling  moment  (L) 
derivatives  differ  substantially  from  the  results  of  case  RK.1-1  and  have  much  larger  Cramer- 
Rao  bounds.  In  addition,  the  calculated  p (unmatched)  does  not  agree  very  well  with  the  measured 
record  (Fig.  5).  The  less  than  good  agreement  even  in  the  matched  records  suggests  an  in- 
adequacy in  the  mathematical  model.  Close  examination  of,  say,  the  sideslip  record  in  Figures 
4 or  S indicates  a possible  non-linearity  with  similar  to  that  noted  in  Reference  4 when  analysing 
the  longitudinal  response.  The  sideslip  angle,  J3,  in  the  present  case  covers  a relatively  wide  range, 
i.e.  approximately  —5°  < < 4°.  The  wind  tunnel  data  of  Reference  8 does  indicate  significant 

changes  in,  for  example,  Nf  and  Lf,  over  this  range  of  P which  would  lead  to  changes  in  the  air- 
craft response  frequency.  In  Reference  4 it  was  possible  to  model  the  non-linearity  of  the 
Pitching  Moment  curve  and  thus  improve  the  match.  In  the  present  case,  a larger  number  of 
parameters  appears  to  be  involved,  and  an  alternative  course  is  adopted.  This  is  to  neglect  that 
initial  part  of  the  record  where  the  sideslip  angle  is  large  and  to  analyse  a later  portion  where 
the  sideslip,  /3,  range  is  smaller  and  non-linearities  less  apparent.  In  Figures  6 and  7 the  first 
1 -85  seconds  of  record  has  been  removed  so  that  the  /3  range  now  varies  over  approximately 
—2-3°  < P <3-7°  with  the  upper  value  only  being  attained  briefly.  In  doing  this  the  rudder 
input  no  longer  appears  so  that  no  rudder  derivatives  can  be  extracted.  In  addition,  the  extrapo- 
lated part  of  the  roll  rate,  p record  (Fig.  3fl)  has  also  been  removed.  In  Figure  6 the  sideslip, 
yaw  rate  and  lateral  acceleration  have  been  matched  but  the  roll  rate  has  been  excluded  while  in 
Figure  7 all  four  records  have  been  included  in  the  analysis. 

The  results  in  columns  RKl-3  and  RKl-4  of  Table  4 together  with  Figures  6 and  7 show  a 
considerable  improvement  in  the  fits  obtained.  Figure  6 shows  almost  perfect  matching  of  the 
sideslip,  yaw  rate  and  lateral  acceleration  records  but  the  agreement  between  calculated  and 
measured  roll  rate  is  still  not  acceptable,  nor  are  the  rolling  moment  derivatives  well  identified. 
In  order  to  improve  these  areas  it  seems  that  the  roll  rate  record  needs  to  be  included  in  the 
matching  process.  When  this  is  done,  the  rolling  moment  derivatives  are  much  better  identified 
as  shown  in  Table  4 (column  RKI-4)  and  the  roll  rate  fit  is  also  considerably  improved  (Fig.  7) 
although  some  room  appears  to  remain  for  further  improvement.  At  the  same  time  the  sideslip 
and  yaw  rate  fits  have  also  improved  marginally  but  the  lateral  acceleration  fit  has  deteriorated 
slightly.  The  derivative  values  identified  do  not  differ  all  that  much  from  those  obtained  in  case 
RKl-1  even  though  the  fits  are  much  better.  This  points  to  the  sensitivity  of  the  time  histories 
to  quite  small  changes  in  some  parameter  values. 

The  results  from  flight  RK2  at  Af  = 0-72  appear  in  Figures  8-10  and  in  Table  5.  The 
rudder  input  used  is  the  smoothed  curve  shown  in  Figure  2b. 


TABLE  5 

Results  from  Flight  RK2 


Item 

A priori 

RK2-1 

RK2-2 

RK2-3 

YglmVe 

-0*292  ±0  058 

-0*292  1 0 054 

-0*308  1 0*055 

-0*290  1 0*055 

YgJmVr 

-0*040  ± 0 008 

-0*043  1 0*008 

-0*042  1 0*008 

— 

yglmV, 

-0*226±  0 045 

-0*1901  0*024 

-0*19410*026 

-0*175  1 0*029 

ygrImVe 

-0*021  ±0*0041 

-0*01810*004 

-0*01910*004 

— 

Lgih 

-25*82±5*2 

-26*4312*46 

-22*3314*93 

-23*5712*55 

Lplh 

-1*78±0*36 

-2*2710*28 

-1*9010*36 

-2*2910*32 

Lrjlx 

0*684  ±0*14 

0*74110*14 

0*690  1 0*14 

0*74110*14 

LgJIx 

-7*57±1*51 

-8*0011*03 

-7*5111*49 

— 

NgUz 

17-44±3*49 

12*5610*32 

12*6010*46 

12*1110*53 

Nplh 

0*01261  0*0025 

0*0126  1 0*0025 

0*0126  1 0*0025 

0*0126  10*0025 

Nrih 

-0*56310*11 

-0*61010*084 

-0*639  1 0*094 

-0*538  1 0*089 

NJlz 

5*7911*16 

5*0410*50 

5*1710*60 

— 

ai 

1*079  1 0*22 

0*9210*12 

0*92±0*12 

0*8610*14 

sHP) 

— 

0*51 

0*33 

0 079 

sHp) 

— 

1*45(0*24) 

— 

0*62(0*10) 

sKr) 

— 

0*42 

0*38 

0 21 

sHg^ylVe) 

0*69 

0*57 

0*25 

The  results  mirror  closely  those  obtained  from  flight  RK 1 and  all  the  remarks  made  previously 
apply  here  also,  although  with  slightly  less  emphasis.  Figure  8 and  column  RK2-1  of  Table  S | 

present  the  results  of  the  full  match  of  all  measured  variables  using  the  model  of  Section  2.3  and  j 

corresponds  to  Figure  4 of  flight  RKl.  A relatively  good  fit  all  round  is  obtained  but  with  some  f 

residual  mismatch  obvious.  The  mean  square  fit  error  for  p in  Table  S is  based  on  the  weighting  | 

element  3906  in  D\  while  the  values  in  brackets  correspond  to  a D\.  of  6S7  (see  comment  following  ^ 

Table  3).  In  Figure  9 the p record  has  not  been  matched  and  this  results  in  some  improvement  in  the 
mean  square  match  error  of  the  other  three  variables  (see  column  RK2-2  of  Table  S)  but  a loss  of 
rollinformationaspreviouslyincaseRKI-2.  Finally,  in  Figure  10  all  the  records  are  once  again 
matched  but  the  first  1 ■ 6 seconds  of  record  removed,  including  the  rudder  input  history  and  the 
roll  rate  extrapolation  of  Figure  36.  The  time  history  fits  are  now  all  good  with  the  corresponding 
parameters  shown  in  column  RK2-3  of  Table  S.  The  values  of  the  extracted  parameters  are 
generally  not  too  different  from  those  of  case  RK2-1  obtained  using  the  full  time  histories. 


4.2  Aerodynamic  Derivative  Estimates 

On  examination  of  the  Cramer-Rao  bounds,  ocrt,  of  the  identified  aerodynamic  derivatives 
in  Tables  4 and  5 it  is  possible  to  classify  the  derivatives  into  three  groups.  In  the  first  group  fall 
those  derivatives  for  which  ocr  shows  no  change  from  the  assigned  a priori  weight,  thus  implying 
little  or  no  information  regarding  these  derivatives  in  the  test  data.  The  derivatives  in  this  group 
are  Lt,  Np,  and  On  the  other  hand,  the  group  of  derivatives  showing  a substantial  change 
in  OCR  are  Nf,  Nt„  Lg,  and  Lg,.  A good  deal  of  confidence  can  be  placed  in  the  identified  values 
of  these  derivatives.  An  intermediate  group  consists  of  Yg,  Lp  and  Nr.  For  these  derivatives, 
some  information  does  exist  in  the  data,  leading  to  some  reduction  of  ocr  relative  to  the  a priori 
weights,  but  they  cannot  be  regarded  as  strongly  identified  as  the  previous  group.  This  classi- 
fication is  in  accord  with  the  groupings  identified  in  Reference  3.  The  point  should  be  stressed 
that  the  above  comments  only  apply  when  all  four  time  histories  are  matched  simultaneously. 
If  the  roll  rate  record  is  neglected  then  the  loss  of  information  severely  reduces  confidence  in  the 
rolling  moment  derivatives. 

Looking  in  more  detail  at  the  results  of  flight  RKl  at  3/  = 0*96  (Table  4),  in  particular 
cases  RKl-1  and  RKl-4  which  match  all  four  records,  it  can  be  seen  that  of  the  strongly  identified 


derivatives,  N^jh,  N,Jlt,  Lflh  and  /-#,//,  only  the  rudder  derivatives  N,JIi  and  L,,llx  are  in 
reasonable  agreement  with  the  a priori  values,  which  are  based  largely  on  wind  tunnel  tests. 
Both  of  the  sideslip  derivatives  Nfjlt  and  £.///<  are  around  20%  less  in  magnitude  than  the  a 
priori  value.  Errors  in  the  moments  of  inertia  /<  or  lx  would  affect  the  rudder  and  the  sideslip 
derivatives  in  a similar  manner  and  hence  do  not  appear  to  explain  the  20%  difference.  A 
similar  comment  applies  to  corrections  for  centre  of  gravity  shift.  Hence  it  appears  that  the 
differences  in  the  sideslip  derivatives  reflect  genuine  aerodynamic  differences  between  flight  and 
wind  tunnel  associated  perhaps  with  the  rapidly  changing  transonic  flowfield  and  perhaps 
coupled  also  with  non-linear  characte  istics  with  sideslip.  The  lower  values  identified  here  should 
by  used  in  any  mathematical  model. 

The  second,  weaker  group  of  derivatives,  Lpjlx  and  Nrlh  also  show  differences 

from  the  a priori  values.  In  all  cases  the  flight  values  are  less  in  magnitude  than  the  corresponding 
a priori  value,  the  difference  varying  from  10  % to  as  much  as  25  %.  As  a result,  the  damping  of  the 
Dutch  Roll  mode,  which  is  proportional  to  the  sum  of  these  three  derivatives,  is  also  well  below 
the  a priori  expectation. 

Of  the  other  parameters  listed  in  Table  4,  the  extracted  values  of  y^ImVe  and  ai  can  be  taken 
with  a certain  amount  of  confidence  because  of  the  decreased  Cramer-Rao  bounds  relative  to 
the  a priori  weights.  The  low  (in  magnitude)  value  of  yflmVe  is  to  some  extent  a reflection  of  the 
extracted  values  of  and  discussed  above,  as  can  be  seen  from  the  definition  of  y^  in 
Equation  (22).  The  value  of  ai  on  the  other  hand,  is  largely  a reflection  of  the  sideslip  vane 
calibration  factor,  G,  as  defined  in  Equation  (19).  Using  this  definition,  the  implied  value  of  G 
is  0-78  for  case  RKI-1  and  0-74  for  case  RKl-4.  For  agreement  with  wind  tunnel  calibrations, 
C should  be  unity.  This  large  reduction  in  C compares  with  the  value  of  0-84  noted  in  Reference 
4 for  the  incidence  vane  at  A/  = 0-96.  While  discussing  the  vane  calibrations,  the  M = 012 
results  (case  RK2)  should  also  be  considered.  The  values  of  ai  from  Table  S lead  to  values 
for  C of  0-86  and  0-80  for  the  two  cases  (RK2-I  and  RK2-3  respectively)  in  which  all  time 
responses  were  matched.  The  equivalent  value  for  the  incidence  vane  at  Af  = 0-71  from 
Reference  4 was  G = 1-01.  Thus  the  present  results  suggest  that  differences  between  tunnel  and 
flight  are  more  pronounced  for  sideslip  measurements  than  for  incidence  measurements  at  both 
subsonic  and  transonic  speeds.  More  work  would  be  required  to  resolve  these  apparent 
differences. 

Turning  now  to  a more  detailed  examination  of  the  identified  aerodynamic  derivatives 
from  flight  RK2  (Table  5)  at  A/  = 0-72,  it  may  be  noted  that  at  this  Mach  number  transonic 
non-linearities  should  not  yet  be  apparent  so  that  a reasonable  agreement  between  experiment 
and  theoretical  estimates  may  be  expected.  Table  6 makes  a comparison  of  the  lateral  derivatives 
(excluding  the  rudder  derivatives)  obtained  from  several  sources. 


TABLE  6 

Comparison  of  M — 0-72  Results 


Derivative 

A priori 

Identified 

Datcom 
(Reference  9) 

Reference  10 

YfImVe 

-0-292 

-0  291 

-0-22 

-0-28 

Lfllx 

-25-82 

-26-43  (-23  -57) 

-33  41 

-24-86 

Lpl/x 

-1-78 

-2-28 

-2-03 

-2-70 

Lrllx 

0-684 

0-741 

0-338 

0-811 

NfHz 

17-44 

12-34 

12-75 

10-83 

NplU 

0-0126 

0-0126 

0-045 

— 

Nrllz 

-0-563 

-0-610  (-0-538) 

-0-235 

-1-12 

The  first  column  in  Table  6 lists  the  a priori  values,  based  mainly  on  wind  tunnel  tests,  for 
case  RK2  at  A/  = 0-72,  while  the  second  column  summarises  the  identified  results  as  in  Table  5. 
These  results  are  the  means  of  columns  RK2-1  and  RK2-3  except  for  those  cases  where  the  two 

13 


I 


columns  differ  by  more  than  10°^  when  both  results  are  quoted.  Column  3 of  Table  6 lists  I 

Datcom  estimates  obtained  from  Reference  9 and  referred  to  flight  dX  M = 0-70  with  lift  co- 
efficient and  c.g.  position  corresponding  to  flight  RK2.  Finally,  column  4 gives  a further  set  of 
results,  appropriate  to  the  Mirage  III.  derived  from  Reference  10. 

There  appears  to  be  reasonable  agreement  between  the  a priori  and  identified  values  of  all 
derivatives  except  for  Lp,  N^,  and  St,  (Table  5).  As  with  flight  RKl  at  A/  = 0-96,  A)//,  is  less 
in  magnitude  than  the  a priori  value  but  in  the  present  case,  unlike  RKl,  Lfjlx  is  in  good  agree- 
ment with  the  a priori  value.  At  the  same  time  the  Datcom  estimate  for  Lf  is  approximately 
30 °a  high  while  the  estimate  from  Reference  10  is  in  good  agreement  (Table  6).  The  present 
value  of  Stilt  is  29%  below  the  a pr/oW  value  while  A,,//,  is  13%  below  (Table  5).  If  it  is  assumed 
that  the  1 3 % difference  in  St,llt  can  be  accounted  for  by  an  error  in  the  moment  of  inertia,  /<, 
this  still  leaves  A^//,  about  16%  in  error.  If,  however,  /*  is  assumed  to  be  13%  in  error  then  the 
theoretical  estimates  for  A^//;  in  Table  6 would  also  be  reduced  by  13%  thus  spoiling  the  rela- 
tively good  agreement  between  flight  and  theoretical  estimates  for  A^.  In  any  case,  there  remains 
a significant  difference  between  flight  and  a priori  (wind  tunnel)  values  for  SflU  at  A/  = 0-72. 

The  other  derivative  of  interst,  Lp,  suggests  a roll  damping  in  flight  almost  30%  greater  than  the 
a priori  estimate,  opposite  to  the  trend  found  at  A/  = 0-96.  The  agreement  with  the  theoretical 
estimates  is  better  although  the  two  theoretical  estimates  differ  by  over  30%.  ^ 

Looking  at  YflmVe,  the  agreement  between  the  various  sources  can  be  regarded  as  reason- 
able but  with  Sr/lt  the  theoretical  estimates  differ  widely  with  Reference  10  giving  a value 
about  twice  the  measured  value  and  Datcom  a value  less  than  half  that.  This  would  lead  to  quite 
large  differences  in  Dutch  Roll  damping.  Finally,  as  may  perhaps  be  expected,  there  are  also 
considerable  differences  between  theoretical  and  measured  values  of  the  weak  cross  derivatives, 

Lr  and  Sp. 


4.3  Fin  Loads  Analysis 

As  a further  example  of  the  use  of  the  modified  Newton-Raphson  procedure,  the  method 
has  been  applied  to  the  estimation  of  fin  load  parameters  from  flight  records  of  Shear,  S,  and 
Tension,  T,  measurements  in  the  fin.  These  records  were  taken  concurrently  with  the  records 
of  sideslip,  roll  rate,  yaw  rate  and  lateral  acceleration  in  flights  RKl  and  RK2  discussed  in  the 
preceding  sections.  The  mathematical  model  used  is  given  in  Section  2.4  and  the  two-part 
procedure  outlined  there  was  used  to  estimate  the  parameters  of  Equation  (29).  The  first  pait, 
involving  analysis  of  the  aircraft  motion  has  been  treated  in  the  preceding  sections  and  the 
identified  parameters  summarised  in  Tables  4 and  5.  Those  parameters  are  now  held  fixed  and 
the  corresponding  Shear  and  Tension  records  analysed  separately.  The  results  are  summarised  in 
Table  7 and  in  Figures  11-14. 


TABLE  7 

Results  of  Fin  Loads  Analysis 


Item 

RKl-l 

RKI-4 

RK2-1 

RK2-3 

Sf 

41-8±9-4 

43-6±lll 

38-3±48-l 

19-81320 

Sp 

-2-0±2-5 

-2-2±3-0 

-6-2±14-9 

-12-4111-4 

Sr 

-7-5±5-9 

-15-3±6-5 

-13-6±24-8 

-20-91171 

S,r 

16-5±5-7 

— 

19-6±5-8 

— 

T, 

88-9±16-9 

95-9±19-9 

150-9±92-4 

91-7161-3 

Tp 

0-27±4-4 

l-06±5-4 

13-6±28-7 

-4-0121-9 

Tr 

-5-2±10-6 

-17-6±ll-6 

16-1147-6 

-10-5132-7 

T,, 

l5-8±10-2 

— 

21-9111-2 

— 

sHS) 

1-34 

0-17 

0-85 

0-47 

sHT) 

0-92 

0-066 

0-54 

0-15 

k 


Columns  RKl-1  and  RKI-4  in  Table  7 correspond  to  the  same  columns  in  Table  4 and  are 
for  flight  at  A/  = 0-96.  The  equivalent  figures  are  Figures  II  and  12  which  should  be  taken  in 
conjunction  with  Figures  4 and  7.  The  fits  in  Figure  1 1 display  the  same  features,  in  particular 
the  apparent  phase  shifts,  as  in  Figure  4 while  Figure  12,  where  the  first  1-85  seconds  of  record 
has  been  removed,  shows  an  excellent  fit  as  did  Figure  7.  Table  7 shows  that,  during  this 
manoeuvre,  the  parameters  most  strongly  identified,  judging  by  the  relative  magnitudes  of  the 
Cramer-Rao  bound,  ocr,  are  the  sideslip  derivatives  and  the  rudder  derivatives  (RK1>1  only). 
There  is  reasonable  agreement  for  the  sideslip  derivatives  between  flight  RKl-1  and  RKM 
but  the  values  of  the  p and  /-derivatives  do  vary  considerably.  It  should  be  remembered  that  no 
I a priori  estimates  have  been  assigned  to  any  of  the  fin  load  parameters. 

• From  the  results  of  RKI-I  and  RKI-4  in  Table  7 it  is  possible  to  estimate  which  derivatives, 

; if  any,  make  the  greatest  contribution  to  fin  loads.  This  can  be  done  by  multiplying  each  derivative 

by  the  maximum  value  of  the  corresponding  independent  variable.  For  example,  the  maximum 
I shear  due  to  roll  rate  would  be  Sp  . pmax-,  etc.  When  this  is  done  it  becomes  apparent  that,  at 

A/  = 0-96  at  least,  the  dominant  contribution  to  both  Shear  and  Tension  comes  from  side- 
slip with  rudder  second  in  importance.  For  Shear  there  is  also  a substantial  contribution  from 
yaw  rate,  r,  but  this  is  less  apparent  with  Tension.  In  both  cases  the  roll  rate  effect  is  small. 

I For  M = 0-72  the  results  are  presented  in  Table  7 under  columns  RK2-1  and  RK2-3  which 

correspond  to  columns  RK2-1  and  RK2-3  of  Table  5.  The  equivalent  figures  are  Figures  13  and 
14  to  be  taken  in  conjunction  with  Figures  8 and  10.  As  previously,  an  improved  fit  is  obtained 
when  the  initial  part  of  the  records,  the  first  1 - 6 seconds  in  this  case,  is  neglected  (compare  Figure 
' 14  with  Figure  13).  From  Table  7 it  appears  that  most  of  the  derivatives  are  weakly  identified, 

with  the  exception  of  the  rudder  derivatives  and  possibly  the  Tension  derivative  Tf.  If,  as 
I before,  the  contributions  to  actual  fin  loads  are  calculated,  the  results  indicate  that  Tension 

>1  forces  are  once  again  dominated  by  sideslip,  /3,  with  roll  rate,  yaw  rate  and  rudder  making  con- 

j siderably  smaller  contributions.  With  Shear,  on  the  other  hand,  no  dominant  contributor  appears 

I and  it  seems  that  all  shear  derivatives  must  be  treated  equally  in  this  case. 

' The  results  obtained  here  must  be  regarded  as  tentative  and  further  work  needs  to  be  done 

i ; (for  example  to  establish  some  sort  of  a priori  structure)  to  increase  confidence  in  them.  Never- 

( ' . theless,  the  analysis  provides  an  example  of  the  sort  of  problem  which  can  be  treated  by  the 

Modified  Newton-Raphson  parameter  estimation  technique. 


5.  CONCLUDING  REMARKS 

A modified  Newton-Raphson  parameter  estimation  procedure  has  been  applied  to  the 
analysis  of  lateral  flight  data  from  a 60°  delta  wing  aircraft  at  two  Mach  numbers,  M = 012 
and  M = 0 -96.  The  response  of  the  aircraft  over  a period  of  5 seconds  {u>nT  = 17  -5)  to  a rudder 
pulse  was  successfully  analysed  to  obtain  all  the  relevant  lateral  parameters.  A data  sampling 
rate  of  60  per  second  (<un^t  = 0 -06)  was  found  to  be  ample  and  no  problems  of  numerical 
convergence  were  encountered.  Records  of  rudder  input,  sideslip,  roll  rate,  yaw  rate  and  lateral 
acceleration  were  used  in  the  analysis  of  aircraft  motion.  Neglect  of  roll  rate  record  was  found  to 
degrade  the  identified  results,  especially  the  rolling  moment  derivatives.  It  was  shown  how 
possible  non-linear  effects  with  sideslip  can  be  avoided  by  basing  the  identification  on  those  parts 


of  the  measured  record  where  sideslip  angle  was  small. 

The  mathematical  model  used  included,  as  an  integral  part,  corrections  to  lateral  acceleration 
measurement  for  lateral  accelerometer  position  error  and  also  included  corrections  to  sideslip 
measurements  due  to  vane  position  error,  vane  dynamics  and  possible  calibration  errors.  The 
results  pointed  to  apparent  differences,  of  considerable  magnitude,  in  sideslip  vane  calibrations 
between  flight  and  wind  tunnel,  which  would  require  further  investigation  to  resolve. 

The  extracted  values  of  the  lateral  derivatives  obtained  in  the  present  study  produced  a 
number  of  differences  from  the  expected  values.  In  particular,  the  yawing  moment  due  to  side- 
slip, Ng,  was  found  to  be  significantly  lower  than  expected  at  both  M = 0-72  and  M = 0-96 
while  rolling  moment  due  to  sideslip,  Lg,  was  low  only  at  M = 0-96.  On  the  other  hand  roll 
damping,  Lp,  was  lower  than  expected  at  Af  = 0-96  but  considerably  higher  at  M = 0-72. 
Each  of  these  effects  is  considered  to  be  largely  due  to  genuine  aerodynamic  differences  as 
opposed  to  errors  in  estimates  of  moments  of  inertia  or  errors  in  centre  of  gravity  position. 


i 

I 


IS 


Some  comparisons  with  theoretical  estimates  was  made  at  A/  - 0*72  with  generally  quite  good 
agreement  apart  from  one  or  two  exceptions. 

A further  application  of  the  method  to  the  analysis  of  fin  loads,  while  producing  only 
tentative  conclusions,  demonstrated  the  possible  adaptation  of  the  identification  procedure  to  a 
variety  of  problems.  In  general  the  modified  Newton-Raphson  technique  has  proved  to  be  an 
effective  and  useful  tool  for  aircraft  flight  test  analysis  with  particular  relevance  to  aircraft 
mathematical  modelling  activities. 
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Coefficients  in  Equation  (19) 

Matrix  of  stability  parameters.  Equation  (2) 

Measurement  bias  vector 

Equation  error  bias  term.  Equation  (24) 

Coefficient  in  Equation  (19) 

Matrix  of  control  parameters.  Equation  (2) 

Unknown  parameter  vector 

Vector  of  a priori  parameter  values 

Weighting  matrix  for  response  variables.  Equation  (1) 

Weighting  matrix  for  parameter  estimates  Equation  (I) 

Matrix  of  stability  parameters.  Equation  (3) 

Gravitational  acceleration 

Matrix  of  control  parameters.  Equation  (3)  or  gain  of  vane.  Equation  (IS) 
Height  above  sea  level,  m 
Moment  of  inertia  in  roll,  kg/m* 

Moment  of  inertia  in  yaw,  kg/m* 

Product  of  inertia,  kg/m* 

Cost  functional  (fit  error) 

Constant  of  proportionality.  Equation  (16) 

Rolling  moment 

Equation  error  bias  term.  Equation  (11) 

Mass,  kg 

Mach  number 

Measurement  noise  vector 

Number  of  data  samples  or  yawing  moment 

Equation  error  bias  term.  Equation  (1 1) 

Lateral  acceleration  in  g units 
Roll  rate,  rad/s 
Pitch  rate,  rad/s 
Yaw  rate,  rad/s 

Reference  area,  m*.  or  Shear  force  in  fin 
Bias  term.  Equation  (28) 


A- 


5*  Mean  square  weighted  error 

t Time,  sec 

T Record  length,  sec,  or  Tension  force  in  fin 

To  Bias  term,  Equation  (28) 

■ Vector  of  control  input  variables 

V Velocity  increment  in  /direction,  m/t 

V Resultant  airspeed,  m/s 

w Velocity  increment  in  z-direction,  m/s 

X Body  axis  co-ordinate  in  forward  direction 

X Vector  of  state  variables 

y Body  axis  co-ordinate  in  lateral  direction 

y Vector  of  response  variables.  Equation  (3) 

V Force  in  >-direction,  N 

Yo  Equation  error  bias  term.  Equation  (I  I) 

Yfy  yp>  Yr,  Yir  Modified  sideforce  derivatives.  Equation  (22) 
Equation  bias  term.  Equation  (24) 

Body  axis  co-ordinate  in  downward  direction 
Measured  response  vector.  Equation  (4) 

Angie  of  incidence  increment,  rad 
Sideslip  angle,  rad 

Inertial  correction  to  p.  Equation  (16) 

Roll  and  yaw  rate  correction  to  fi,  Equation  (14) 
Increment 
Aileron  angle,  rad 
Elevator  angle  increment,  rad 
Rudder  angle,  rad 
Pitch  attitude  increment,  rad 
Air  density,  kg/m* 

Standard  deviation 
Cramer-Rao  Bound 
Roll  attitude,  rad 

Natural  undamped  frequency,  rad/s 

Subscripts 

a Lateral  accelerometer 

e Trim  or  equilibrium  state 

/ Time  index 

fn  Measured  value 


Derivative  w.r.t.  p 


Yo 

z 

z 

a 

P 

P4 
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A 

8« 

S« 

8r 

e 

p 

a 

OCR 
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Derivative  w.r.t.  r 


Sideslip  vane 
Derivative  w.r.t.  /S 
Derivative  w.r.t.  ^ 


Derivative  w.r.t.  fi* 


Derivative  w.r.t.  8, 
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